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Abstract 

The scaling of decoherence rates with qubit number N is studied for 
a simple model of a quantum computer in the situation where N is large. 
The two state qubits are localised around well-separated positions via 
trapping potentials and vibrational centre of mass motion of the qubits 
occurs. Coherent one and two qubit gating processes are controlled by 
external classical fields and facilitated by a cavity mode ancilla. Deco- 
herence due to qubit coupling to a bath of spontaneous modes, cavity 
decay modes and to the vibrational modes is treated. A non-Markovian 
treatment of the short time behaviour of the fidelity is presented, and 
expressions for the characteristic decoherence time scales obtained for the 
case where the qubit/cavity mode ancilla is in a pure state and the baths 
are in thermal states. Specific results are given for the case where the 
cavity mode is in the vacuum state and gating processes are absent and 
the qubits are in (a) the Hadamard state (b) the GHZ state. 



1 C Introduction 



Quantum computers are expected to have a key advantage over classical comput- 
ers in terms of computational complexity, an advantage based on the idealised 
features of parallelism, entanglement and unitary quantum computation pro- 
cesses. For certain algorithms 0111, the number of computational steps needed 
using quantum computers should increase much more slowly with input number 
size than for classical computers, enabling computations that were infeasible us- 
ing classical computers to be carried out. In the case of the quantum search 
algorithm for example, the number of steps increases with the square root of 
the number of items, rather than linearly as in classical searching However, 
the physical system whose quantum states define the N qubit system and the 
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quantum devices involved in the gating processes both interact with the envi- 
ronment. Such interactions change the density operator describing the quantum 
computer state from a pure to a mixed state, with coherences between states 
evolved from different input states being partially or completely destroyed. This 
process of decoherence is the enemy of quantum computation, though methods 
such as quantum error correction [310], decoherence free sub-spaces [51ini[7||S] 
and dynamical suppression of decoherence [51 110| could be used to minimize 
its effects. Utilizing the computational complexity advantage of quantum com- 
puters requires suitably large qubit numbers (ca 10^ qubits may be needed for 
searching or factoring algorithms with error correction for input numbers 
of practical importance), so the scaling with the number of qubits of the de- 
coherence rates in comparison to the coherent processing rates is important in 
determining the size limits of useful quantum computers . 

This paper extends previous work El El El El Ej treating scaling 
effects. Decoherence effects are studied for a simple model of an N qubit quan- 
tum computer, involving TV two state qubit systems (see figure 1). The qubits 
are localised around well-separated positions via trapping potentials, and the 
centre of mass (CM) vibrational motions of the qubits are treated. Coherent 
one and two qubit gating processes are controlled by time dependent localised 
classical electromagnetic (EM) fields and magnetic fields that address specific 
qubits. The two qubit gating processes are facilitated by a cavity mode ancilla, 
which permits state interchange between qubits. The magnetic fields are used 
to bring specific qubits into resonance with the classical EM fields or the cavity 
mode. The two state qubits are coupled to a bath of EM field spontaneous 
emission (SE) modes, and the cavity mode is coupled to a bath of cavity decay 
modes. For large N the numerous vibrational modes of the qubits also act as 
a reservoir, coupled to the qubits, the cavity mode and the SE modes. The 
system-environment coupling interactions include electric dipole coupling of the 
qubits to the SE modes, quasi-mode coupling of the cavity mode to cavity decay 
EM field modes and Lamb-Dicke coupling of qubits to CM vibrational modes 
and each of: (a) the SE modes (b) the cavity mode (c) the classical EM gating 
fields (d) the classical magnetic gating field. The system-environment coupling 
interactions considered include both az phase destroying and ax population 
destroying terms, though the latter are more important in this model. The 
model is similar to those of the ion-trap pi!, neutral atom 22 and cavity QED 
|23ElEj varieties. 

The object is to study decoherence effects in quantum computers for the 
situation where the number of qubits N becomes large. Decoherence effects will 
be specified by the fidelity, which measures how close is the actual behaviour of 
the density operator for the qubits system (including the ancilla) to its idealised 
behaviour due to coherent gating processes only. The initial quantum state for 
the quantum computer will be assumed to be pure and the reservoir states will 
be assumed to be thermal. Previous work (HI El El El El El E] indicates 
that the decoherence time tD can decrease inversely with N (qubits decohering 
independently) or N'^ (qubits decohering collectively), so it can become very 
short for the case where N is large and may become comparable to the reservoir 
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correlation time tc- In the case of atomic spontaneous emission relaxation, tc 
~ 10"^'' s, and ~ 10~^s for one qubit could become to ~ 10~^^s for 10^ 
qubits in the case of collective decoherence. However, one aspect of the present 
study is whether the qubits decohere independently or not, so we do not arbi- 
trarily assume either case. As maintaining coherence between all states of the 
quantum computer is important for its operation, we are interested in the be- 
haviour of all density matrix elements, and not just those for a smaller density 
matrix describing a single qubit, even though the latter may have long deco- 
herence times and satisfy Markovian equations. Treating the density operator 
for the full N qubit system via the standard Born-Markoff master equation (see 
for example, |26[ I27|) requires that all the interaction picture density matrix 
elements do not change significantly during the reservoir correlation times. As 
the latter are just the time scales over which two-time correlation functions of 
reservoir operators involved in the system-reservoir interaction decay to zero, 
we can see that the reservoir correlation times are independent of the number 
of qubits. Since the general decrease in decoherence time with increasing qubit 
numbers indicates that some density matrix elements must change over shorter 
and shorter time scales, then using Markovian master equations for studying 
decoherence therefore becomes questionable if the number of qubits becomes 
large enough In the present paper we do not assume that Markovian behaviour 
will necessarily occur, and use methods that do not rely on the standard Born- 
Markoff master equation. If the behaviour does turn out to be Markovian then 
our approach is still correct. However as we will see, in certain cases the be- 
haviour is not consistent with the Markoff approximation, so the necessity for 
our choice of non-Markovian methods is justified a posteriori. Several authors 
j28l 1291 1801 m] have also emphasised the importance of studying the short time 
regime for macroscopic and mesoscopic systems (such as quantum computers 
with large qubit numbers) using non-Markovian methods. One such approach 
PnilSni involves using eigenstates of the system operators involved in the system- 
reservoir interaction (pointer basis), another treatment [221 develops short-time 
expressions for the linear entropy (or impurity). Here the approach used is 
similar to that in |33| and uses the Liouville-von Neumann equations to ob- 
tain short-time expansions for the fidelity. This leads to expressions for the 
characteristic decoherence time scales for the general case with gating EM and 
magnetic fields present, and at non-zero temperature, and with the qubit and 
ancilla system in an arbitrary pure state. The treatment enables situations 
where independent or collective decoherence occur to be distinguished. In this 
initial paper, the decoherence time scales will be evaluated for specific initial 
quantum computer states such as: (a) the equal superposition of all qubit states 
(Hadamard state) and (b) the generalised GHZ states, with the cavity mode an- 
cilla in the vacuum state and gating processes ignored. The effects of gating 
processes and the relationship of the results to those obtainable for longer times 
using Markovian master equation methods will be treated later. 

In Section 2 the general theory used will be outlined. Results for the deco- 
herence time scales are presented in Section 3 and Section 4 contains a summary 
of the paper. 
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2 Theory 



2.1 Hamiltonian 

The total Hamiltonian is written as the sum: 

H = Hs + Hc + Hb + Vs + Vi (1) 

The Hamiltonian for the qubit system and cavity mode ancilla is given in 
terms of qubit atomic spin operators cr^ ^-iid cavity mode annihilation, creation 
operators 6, 5^ by: 

Hs = J2 ^^o4 + hujbb^b (2) 

i 

Here the qubits all have transition frequency loq and Ub is the cavity mode 
frequency. The cavity mode facilitates state transfer between different qubits in 
two qubit gating process. In the case of a large number of charged qubits, the 
CM vibrational modes are too closely spaced to enable the in-phase vibrational 
mode to be used in this role, as in the case of ion-trap quantum computers |21| . 
The array of qubits confined in their trapping potentials is rather analogous to 
a low density solid. Also, because the cavity mode energy quantum is large, this 
ancilla is not subject to thermal effects as are the low frequency CM vibrational 
modes. 

The Hamiltonian for the centre of mass vibrational motion of the qubit 
system is given in the harmonic approximation as: 

ia ija.f3 

= J2 fi'^KA^'j^AK (4) 

K 

In these equations the qubit CM momentum operators are pia and the centre 
of mass of each qubit vibrates with displacements 6ria about well-separated 
positions r^o defined by trapping potentials = x,y,z). The vibrational 

mode K has angular frequency and phonon annihilation, creation operators 
Ak, A\^. The vibrational Hamiltonian He applies to the general case where the 
qubits vibrate collectively, such as in the ion-trap case where the charged qubits 
interact via long-range Coulomb potentials. If the qubits vibrate independently, 
such as for neutral qubits with dipole-dipolc couplings ignored, then the coupling 
matrix is given by V^^ — SijSafsV, and the vibrational modes K are specified 
ia. These all have the same frequency v. 

The Hamiltonian for the bath of spontaneous emission and cavity decay 
modes is of the form: 

iJs = E ^^ko-lak + E ^^^kbih (5) 

k k 
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Here the SE mode k has angular frequency ujk and annihilation, creation oper- 
ators afc, a|,, and the cavity decay mode k has angular frequency ^j. and annihi- 
lation, creation operators 6fc,6^. These modes are described in the quasi- mode 
picture [51] . 

The coherent coupling (or gating) processes for the qubit system are de- 
scribed by the Hamiltonian: 

Vs = Y.^hn,a\+Hc) + Y,\m^l~ ^^oyz 

i i 

+ Y.{hg^a'xb + Hc) (6) 

i 

In order, the terms are the coupling of the qubits to the classical EM field, to 
the magnetic field and to the cavity mode. Each qubit can be addressed by 
localised classical EM fields E and magnetic fields B to facilitate one and two 
qubit gating processes. The two qubit gating process involves the creation and 
annihilation of a single photon in the cavity mode, which acts as an ancilla. 
The quantities fli, AiaSLnd gi are coupling constants. The coupling constants 
Jli, are both spatially localised and time dependent in accordance with the 
sequence of gating processes involved for a particular algorithm. 

The interaction of the qubit system and ancilla with the baths and the qubits 
centre of mass vibrational degrees of freedom is given by: 

+ He) + Y^{hh\ukbk -f Wkbl) + He) 

ik k 

+ Y.^f^'^kK^xak[AK + A^k\ + He) 

ikK 

+ Y^{hp'Ka'xb[AK + A^K\^He) 

iK 

+ Y,{f^Q^K^x[AK + A^k]+ He) 

iK 

+ E ^^KkII - a^z) + (1 + <^'z))[Ak + A^k] (7) 

iK 

The first two terms are the electric dipole coupling of the qubits to SE modes 
and the quasi-mode coupling of the cavity mode to the cavity decay modes. 
These terms involve coupling constants g\,Uk and Wk- The remaining terms 
are Lamb-Dicke couplings with the various fields, allowing for the vibrational 
motion of the qubits CM around their equilibrium positions rjQ. The third 
term allows for coupling to the SE modes, the fourth to the cavity mode, the 
fifth to the classical EM field and the sixth to the magnetic field. These terms 
involve coupling constants "■fci^-jPjj-, "^oif- Since this paper studies 

short time scale effects, the standard rotating wave approximation (RWA) has 
not been made. Other possible interactions, such as those due to Roentgen 
currents j85l Kffij for neutral and charged qubit systems, and those due to ionic 
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currents [3311201 ^'^^ charged qubit systems, were also investigated but found to 
be relatively small in comparison to the interactions listed above. 



2.1.1 Expressions 

Detailed definitions or commutation rules for the qubit, cavity mode, CM vi- 
brational modes and bath modes operators are as follows: 



'X 



(|1)(0| + |0>(1|). aV = -^(|l)(0|-|0)(l|). 
<^z = (|1)(1|-|0)(0|). (8) 
[6,5t] ^ 1 [Ak,AI]=5kl 
[ak,aj] = Ski [bk,bl]^Ski (9) 

The qubits are represented by Pauli spin operators, whereas bosonic annihilation 
and creation operators apply for the cavity mode (quantized in volume Vb), the 
CM vibrational modes, the SE modes and the cavity decay modes (quantized 
in volume V). 

The CM displacements Sria {a — x, y, z) are related to the vibrational nor- 
mal coordinates via a unitary matrix S in the form: 



Sna = J2 ^-^^V^^^^^' + ^A'): (10) 

where the matrix S is determined from the eigenvalue equations 

^ifSj0;K = rnv\S.,a,K- (H) 

These equations also give the vibrational frequencies vk- For the case where 
the qubits vibrate independently Sia-^K — Sia-K and ly — {V/m)^. 
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Explicit expressions for the coupling constants are: 



fli = -i^J-^^^{dio-ec)acexpi{\s.c-rio-uJct) (12) 



2€ohV 

= -[vCVaa-^]/h a = 0, 1 (13) 



9i = 2^JiVb ' ®^P(*'*b ■ '"io) (14) 



9k = -«J^^^(dio •efc)exp(ik-rio) (15) 





1 ^ 


2€ohV\ 


1 2muK 



X exp(ik • Tio) (16) 





1 ^ 




1 2mvK 



X cxp(ik6 • Tjo) (17) 



X ac exp i(kc • Fjo - Wei) (18) 

dRif} V 21711/ K 



<K = -2^i^^^-^Ul7^_S.P;K /h a = 0,1. (19) 



Here dio is the electric dipole matrix element, maa(a = 0. 1) arc the magnetic 
dipole matrix elements for each qubit. The polarisation unit vectors for the 
classical field, the cavity mode, the SE mode or cavity decay mode k are ec, e;,, 
respectively, and kc,kh,k give the wave vectors. The classical EM gating field 
could be thought of as being represented by a multi-mode coherent state, where 
ac is the coherent state amplitude for the mode c. All the modes have essentially 
the same polarisation vector, wave vector and frequency, though the frequency 
and wave vector bandwidths are sufficiently wide to produce the required spatio- 
temporal localisation. The equivalent classical electric field at the equilibrium 

position of the i qubit is given by E = i J^c ^c^c exp i(kc • Fjo — uJct) + cc. 



2.2 Dynamics 

The total density operator W for the full system of qubits, qubits CM motion, 
cavity mode and baths satisfies the Liouville-von Neumann equation: 

dW 

ih— = [H,W] (20) 

Initial conditions are assumed of an uncorrelated state for qubits and ancilla, 
bath and CM vibrational motion given by:. 

W^(0)=Ps(0)pb(0)pc(0) (21) 
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and for the bath and CM in thermal states the average value of Vj is zero. 

TrBcViPBiO)PciO) = TrBcPBiO)PciO)Vi = (22) 

In the present paper the initial state for the qubits and ancilla will be assumed 
in a pure state so that ^5(0) = |V's)('/'sl- 

The dynamics of the quantum computer is described by the reduced density 
operator pg for the qubits and ancilla subsystem defined by: 

Ps - TrscW (23) 

The exact evolution allowing for the coherent coupling and the interaction with 
the bath and the CM vibrational modes is embodied in the time dependence of 
Pg. No assumption is made here that pg satisfies a Markovian master equation, 
though it may do so for time scales which are large compared to bath correlation 
times. 

The coherent evolution of the qubits and ancilla subsystem is described by 
the reduced density operator p^Q, which is chosen to coincide with the exact 
reduced density operator at time zero, and is then allowed to evolve under the 
coherent coupling term Vs only. Hence: 

= [Hs + Vs,Pso] (24) 
Pso(O) = Ps(0) (25) 

Thus the interaction between the qubits/ancilla and the baths and CM motion 
is considered to be switched off here, hence pgQ evolves according to a Liouville- 
von Neumann equation.. 



2.3 Decoherence 

Decoherence effects are specified by the fidelity, defined as: 

F^TrsipsoPs) (26) 

The fidelity specifies how close the actual evolution of qubits and ancilla is to 
the idealised coherent evolution. The time dependence of the fidelity is entirely 
due to the decoherence effects caused by the bath and CM interactions with the 
qubits and cavity mode ancilla. 

As in |33| the short-time behaviour of the fidelity with the qubits in a pure 
state can be expressed as a power series in the time elapsed and explicit expres- 



8 



sions obtained for the first few time constants involved. Thus: 




(27) 



h 







(28) 




TrBc{i's\ViiOf\i's)pBiO)PciO) 

-TrBc{i^s\ViiOMsfPB{0)Pc{0) 
((AViiO)')s)BC 



(29) 
(30) 



The initial value of the fidelity is unity. The times ti,T2, ..specify characteris- 
tic decoherence times for the qubit and ancilla system, their inverses defining 
dccohcrcncc rates. For thermal bath and CM states, only T2 is involved in spec- 
ifying short time decoherence. For the time constant T2, the expression involves 
the average of the square of the fluctuation AV7(0) = Vj(0) — {Vi{0))s of the 
zero time interaction operator Vf(0). The squared fluctuation operator is first 
averaged over the initial qubit and ancilla pure state IV's)) then the result is 
averaged over the bath and CM initial state Pb(0)Pc(0)- 

3 Results 

3.1 General case 

The evaluation of T2 for the general case with gating EM and magnetic fields 
both present, and at non-zero temperature, and with the qubit and ancilla 
system in an arbitrary pure state, leads to expressions for T2 that involve: 

(a) Expectation values for the state \ipg) of one qubit operators uj^, two 
qubit operators ct^ct^ (i ^ j), (a, (3 = X, Z). 

(b) Expectation values of cavity operators b, 6^, bb^ , b^b, b^, 6'''^. 

(c) Expectation values of products of these qubit and cavity operators. 

(d) Sums over bath modes and centre of mass vibrational modes of prod- 
ucts of pairs of the coupling constants 3^, Wfej ^^feX'^'K' '^aif ' ^^'^ products 
being weighted by factors involving the thermally averaged photon and phonon 
quantum numbers nk,fnk, Nk for the SE modes, cavity decay modes and vi- 
brational modes. The thermally averaged quantum numbers are given by the 
Planck function. 

The general case will be studied in further work. 

3.2 Case of spontaneous emission for stationary qubits, 
no cavity mode 

In this case the cavity mode is ignored, the qubits are assumed stationary so that 
CM vibrational modes are excluded (as are all Lamb-Dicke interactions), and the 
only coupling constants included are g^, given by gk expik ■ Vio- The qubit pure 
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state is |0q). This special situation has been studied in where temperature 
effects and criteria for the qubits decohering independently or collectively were 
examined. 

For this case the decoherence time ti is given by: 

^ = E((^^^'^)-('^^)('^x»E5fc5r(2^fe + l) (31) 

= E(^'^xAc7^x)El5fc|'cosk.d,,coth^ (32) 

y fc 

where (..) = ((/<q|..|(/)q), Af2 = J7 — (f2) and d^j = r^o — fjo- This result is the 
same as in P^. 

As in 33 , if Afc is the bandwidth and fc the mean of a normalised Gaussian 
model for the function j g^p coth -^jr^, then the condition for independent deco- 
herence (where only terms with i — j contribute to T2) requires that Akdij ^ 1 
for any pair of qubits {dij — \dij\). Collective decoherence (where terms from all 
pairs I 7^ J of qubits contribute to T2) requires that Akdij ^ 1 and kdij <^ 1. 
Further discussion about independent and collective decoherence is given be- 
low for the case of zero temperature, but based on evaluating the sum over SE 
modes fc for a three dimensional model. 



3.3 Case of no gating processes, zero temperature and 
cavity in vacuum state 

The case where no gating is taking place, the ancilla is back in its original no 
photon state |0)yi and the temperature is zero, provides a simple case of the 
general results, enabling the dependence on the qubit state to be examined. 
Here — |(/)q)|0)a, where the state of the qubit system is |0q). This case is 
of some physical interest, since it corresponds to states produced after idealised 
coherent gating processes have occurred. Too rapid a decoherence for such states 
would be of concern for the implementation of quantum computers. 

3.3.1 Case of correlated qubit states 

ij 

+ Y.^o\a'x)K^,+L (33) 
Y.9\^€+Y.<K<K (34) 

fc kK 

Y.^kp'k (35) 

K 

fc 



2ri 



L = 
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For the correlated qubits case, the decoherence time depends on SE coupling 
constants, non-RWA cavity mode coupling constants and Lamb-Dicke coupling 
constants involving the cavity mode and the SE modes. Note that terms involv- 
ing the same qubits (i = j) and different (i ^ j) qubits are involved. 



In this case {a^a-'y) — for i ^ giving the decoherence time via: 



3.3.2 Case of uncorrelated qubit states 

= Y.^l-\{a'x)?){Gu+K,. 

i 



1 

27| 



(37) 



E(i - i(-x>n(E \9i\" + E Hk\" + E \pk\") 



kK 



K 



+ L 



(38) 



For the uncorrelated qubits case, the decoherence time depends on the same 
coupling constants as for the correlated case, but now only diagonal terms in- 
volving the same qubits are involved. 



3.3.3 Special sub-case of Hadamard or fiducial state 

This particular uncorrelated state results from applying the Hadamard one qubit 
gating process to every qubit, and which produces the important fiducial state 
of an equal superposition of states representing all numbers from to 2^ — 1. 
The state vector for the i qubit is I^q)^ = (|0)i -I- |l)i)/2^/^. In this case the 
decoherence time is given by: 



2ri 



PkPk 

ijK 



ijK k 



(39) 
(40) 



The first term is associated with the characteristic time scale due to the Lamb- 
Dicke term coupling the qubits, cavity mode and vibrational motion. For the 
general case of collective qubit vibrations, this term is only evaluated approx- 
imately, the l/i'K factor appearing in the sum over K being replaced by an 
average value, and thus enabling the unitary properties of the matrix S to be 
employed. Only i — j terms then result. For the case of independent qubit 
vibrations, result H4U|I is exact. The Lamb-Dicke parameter 77 « k • <5r is small 
compared to unity, « (jj^^v^)^^^ l^^iol is the cavity vacuum Rabi frequency. 
The second term is associated with the characteristic time scale due to the 
non-RWA term coupling the cavity mode and cavity decay modes. Ignoring 
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this term, the decoherence time scales as l/\fN for the Hadamard state. For 
Q}) w 50 Mhz, rj ~ 10~^ and thus -qgb ~ 10^s~^, we find that the decoherence 
time scale is about 10~^s for A'' w 10^ qubits. If the quantum computer were to 
become truly macroscopic with N « 10^^, the decoherence time is about lO^^^s. 
Decoherence effects associated with the SE modes do not occur for this state, 
as can also been seen from the expressions for the case of stationary qubits. 

3.3.4 Special sub-case of GHZ state 

This correlated state is of special interest due to its highly entangled nature. 
The state vector for the qubit system is = (|00..0) + |ll..l))/2i/2and the 



ik ikK 

+Y.\^K? +Y.\''k? (41) 

iK k 

NY.\9k? (42) 

fc 

Nr}'Y.\9^\' (43) 

k 

Nrfgl (44) 

The first term in the expression for T2 is associated with the characteristic time 
scale due to the electric dipole term coupling the qubits and SE modes. It can be 
expressed in terms of the RMS vacuum electric field E^tos for the SE modes as 
proportional to the square of the Rabi frequency (Erms ■ dio/ft) associated with 
Erms- The vacuum Rabi frequency for the k SE mode is gk ~ (26^^^)^^'^ |dio|- 
The sum J2k IS'fcP t)e evaluated for SE modes in three dimensions by incor- 
porating a cut-off factor for Ipfep of the form cxp{—LUk/^c), where lOc is cut-off 
frequency. The correlation time for the bath of SE modes is given by the inverse 
of the cut-off frequency and is about 10~^''s. We find that: 

El5fcl'«rwc(c^>o)' (45) 

k 

Here T is the SE decay rate. For F w lO^s'S wq ~ lO^^s''^, ojc « lO^^s'^ we 
find that (^j, l.gfep)^/^ « 10'^'''s^^.Thc second term in the expression for T2 is 
associated with the characteristic time scale due to the Lamb-Dicke term cou- 
pling the qubits, the SE modes and the CM vibrational motion. For the general 
case of collective qubit vibrations, this term is only evaluated approximately. 
For the case of independent qubit vibrations, the result is exact. As the Lamb- 
Dicke parameter is small compared to unity, the second term can be ignored 



decoherence time obtained from: 

1 

ik 

Ek^i' « 

ikK 
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in comparison with the first. The third term in the expression for T2 is associ- 
ated with the characteristic time scale due to the Lamb-Dicke term coupling the 
qubits, the cavity mode and the CM vibrational motion. For the general case 
of collective qubit vibrations, this term is also only evaluated approximately. 
For the case of independent qubit vibrations, the result is exact. It involves the 
square of the Rabi frequency associated with the cavity mode. For gi, « 50 Mhz, 
77 w 10~^ and thus rjgb w 10^s~^, so the third term is also small compared to the 
first. The fourth term in the expression for T2 is associated with the non-RWA 
term coupling the cavity mode and cavity decay modes. Ignoring the last term, 
the decoherence time scales as for the GHZ state. As the first term is 

the most important, we see that the decoherence time scale is about 10~^''s for 
N w 10^ qubits. If the quantum computer became macroscopic with A'' w 10^^ 
the decoherence time is even shorter, about 10~^^s. Evidently the GHZ state 
would not preserve its entanglement for very long. 

3.3.5 Case of correlated qubit states - spontaneous emission deco- 
herence only 

The question of whether the qubits decohere collectively or independently is 
now examined further. For the case of stationary-qubits at zero temperature 
and with no cavity mode ancilla and only the coupling to the SE modes included, 
the decoherence time T2 is given as a special case of the previous expression as: 

^ = ^(Aa^Aa^^) ^ cosk • d,, (46) 

Evaluating the sum \gk\'^ cosk • dij for SE modes in three dimensions and 
incorporating a cut-off factor for Igkl"^ of the form exp(— w^/wc) as before, the 
sum is given by the expression: 

^ l^fep cosk • dij « rujc{ujc/uJofF{ujcTtj, cos^ 6'y;io) (47) 

k 

where r^j — dij/c is the time for light to travel between the i and j qubits, 
Oij;io is the angle between the vectors for the separation of the qubits dij 
and the dipole moment dip. The function F is of order unity and involves 
ratios of polynomials in the quantity uJcTij. For ujcTij 2> 1 this function 
behaves like (wcTy)~^. This is in contrast to the one dimensional case at 
low temperatures treated in , ,33; . where no such decrease for large uJcTij oc- 
curs. For r « 10*s^\ Wo ~ 10^^s~\ ujc ~ 10^ ''s~^ (as above) we find that 
{ruJc{uJc/i^o)^y^^ ~ lO^^s"^. Allowing for conveniently addressing separate 
qubits with optical laser fields, the nearest neighbour distance will be taken as 
dij w lO^^m, hence ujcTij ~ 10"^ ^ 1. Thus even for nearest neighbour qubits, 
{Tujc{ujc/ ojqY Py/'^ w lO^s"^ . This is a factor of 10^ smaller than for the i = j 
terms. For more widely separated qubits the reduction is even greater. The 
greater number of qubits in a shell of qubits between dij , dij -\- Sdij (proportional 
to df ) does not overcome the d"'* factor. Thus it would seem that only the i = j 
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terms will be important in the expression for the dccohcrcncc time T2 . This in- 
dicates that the qubits decohere independently. Including then only the i = j 
terms (which are all equal), the expression for 1/2t2 is proportional to N and 
is non zero for correlated states. Overall, the dccoherence time scales as 
for correlated states. Similar considerations to before again lead to a dccoher- 
ence time scale of about 10~^^s for N w 10^ qubits, whilst for a macroscopic 
size quantum computer with N « 10^^ the dccoherence time is about 10~^^s. 
Remembering that the SE bath correlation time is of order ojj^ 10~^^s, then 
these dccoherence times are not long compared to the correlation time, so non- 
Markovian conditions apply as has been expected. If spontaneous emission is 
present this would seem to be something of a problem for quantum computation 
based on the model treated here, since processing involving correlated states is 
an essential feature. 

3.4 Case of no gating processes, zero temperature, cavity 
in vacuum state but spontaneous emission absent 

The case where no gating is taking place, the ancilla is back in its original no 
photon state |0)^, the temperature is zero, and SE modes arc ignored provides 
a simple case of the general results, enabling the dependence on the qubit state 
to be examined. Here it is assumed that the cavity is of such high quality that 
qubit decay via the SE modes can be ignored, so the only decay is via Lamb- 
Dicke coupling through the cavity mode. Again \ipg) = |0q)|O)a, where the 
state of the qubit system is \(f>Q). This case is of some physical interest, since 
it corresponds to states produced after idealised coherent gating processes have 
occurred, but now in a more favourable cavity situation. Too rapid dccoherence 
of correlated qubit states - even in this case where SE decay is absent - would 
be of concern for the implementation of quantum computers. 

3.4.1 Case of Correlated qubit States 

The previous expressions can be used to obtain the dccoherence time T2 by 
deleting all contributions associated with the SE modes. This gives: 

J2{<^'^<^'x)Kij + L (48) 

ij 

Y.PkP'k L = Y,w,wl (49) 

K k 

3.4.2 Evaluation of Lamb-Dicke Term 

For the general case of collective qubit vibrations, an approximate evaluation of 
the Lamb-Dicke term yields the result: 

i « Nr^^gl (50) 



2r^ 
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Result (I5UII is exact for the case of independent qubit vibrations. This term 
is associated with the characteristic time scale due to the Lamb-Dicke term 
coupling the qubits, cavity mode and vibrational motion. The Lamb-Dicke 
parameter « k • is small compared to unity, gt, « (2^^^)^^'^ \^w\ is the 
cavity vacuum Rabi frequency. The decoherence time scales as 1/ ^/N . For gi, w 
50 Mhz, rj » 10"^ and thus rjgt « 10^s~^, the decoherence time scale is found 
to be about lO^^s for N « 10^ qubits, whilst for a macroscopic sized quantum 
computer with N « 10^^ the decoherence time is about 10~^^s. Although it 
might be thought that these decoherence times are rather long compared to 
bath correlation times, it should be noted that for the CM vibrational modes 
bath, the correlation time is of order where v is a typical CM vibrational 
frequency. As ly ^ 27r.lO® s~^ ,21,,^, the correlation time for the bath is about 
10~^s. Thus these decoherence times are short compared to the correlation 
time, so non-Markovian conditions apply as has been expected. The first result 
is more promising from the point of view of implementing quantum computers 
of a useful size, for example with N w 10"*. In the general case where the 
qubits vibrate collectively, the evaluation is only approximate however and a 
full treatment will involve determining the matrix elements Sia;K relating qubit 
CM displacements to the vibrational normal coordinates in order to examine 
issues such as whether the qubits decohere independently or collectively. 

4 Summary 

Decoherence effects in quantum computers have been studied for the situation 
where the number of qubits N becomes large. A standard model involving N 
two state qubit systems was treated, with localised, well-separated qubits un- 
dergoing vibrational motion in the trapping potentials. Coherent one and two 
qubit gating processes were controlled by time dependent localised classical EM 
and magnetic fields, the two qubit gating processes being facilitated by a cavity 
mode ancilla. The qubits were coupled to bath of spontaneous emission modes, 
the cavity mode was coupled to a bath of cavity decay modes. The numerous 
vibrational qubit modes also behave as a reservoir. Non-RWA couplings and 
effects of qubit vibrational motion were both included. The predominant cou- 
pling of qubits to the baths is amplitude coupling via the a\ Pauli operators. 
Decoherence effects were specified by the fidelity, and non-Markovian expres- 
sions for the short time behaviour of the fidelity were obtained for the general 
case where the qubits and ancilla are in any pure state and the baths and qubit 
vibrational modes are in thermal states. 

Characteristic decoherence time scales were evaluated for specific qubit states 
in the simple situation when no gating processes are occurring and at zero tem- 
perature. The decoherence time is mainly affected by the spontaneous emission 
decay of the qubits. Decoherence times scaling inversely as were found 

for the cases treated. For the uncorrelated Hadamard state a decoherence time 
scale of about lO^^s for N « 10* qubits was found, for a macroscopic N w 10^^ 
it was about 10^ ^^s. For the correlated GHZ state a decoherence time scale of 
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about 10~^''s for N w lO'^ qubits was found, for a macroscopic size quantum 
computer with N » 10^^ it was about lO^^^s. Correlated states generally were 
found to have similar decoherence times to the GHZ state, and the qubits were 
found to decohere independently. 

Characteristic decoherence time scales were also evaluated for specific qubit 
states in a further simple situation when no gating processes are occurring, 
the temperature is zero, and decay via spontaneous emission modes is ignored. 
The decoherence time is now mainly affected by the Lamb-Dicke coupling of the 
qubits to the cavity mode ancilla and to the CM vibrational modes. Decoherence 
times scaling inversely as for the cases treated were found. A decoherence 

time scale of about 10~®s for N « 10'' qubits was found, for a macroscopic size 
quantum computer with N « 10^^ it was about 10~^^s. 

The case where spontaneous emission decay can be disregarded, such as when 
the qubits are all in a very high Q cavity, is far more promising from the point 
of view of implementing quantum computers of useful size, say with A'^ w 10^ 
qubits. An alternative, though more complex model for avoiding spontaneous 
emission effects would be to replace the two internal state qubits by qubits with 
three internal states, such as in a lambda system. In such a model the qubit 
states would be two near degenerate ground states, and the upper state would 
be in effect part of the ancilla system and only coming into play during the 
gating processes. Spontaneous emission would then only be important during 
the gating processes and could be minimized if the latter occurred in short 
enough times. A similar model has been examined in a different context in [7j. 

Finally, although the decoherence times found here are not long compared to 
bath correlation times, and hence non-Markovian conditions apply as expected, 
it would be of interest to consider the behaviour of the fidelity in the Marko- 
vian regime for the present model in order that the short time scale results for 
decoherence effects presented here can be linked up with those for longer time 
scales. 



5 Figure caption 

Figure 1. Model of an A^ qubit quantum computer. Two state qubits are lo- 
calised around well-separated positions via trapping potentials, and undergo 
centre of mass (CM) vibrational motions. Coherent one and two qubit gating 
processes are controlled by time dependent localised classical electromagnetic 
(EM) fields and magnetic fields that address specific qubits. Two qubit gating 
processes are facilitated by a cavity mode ancilla, which permits state inter- 
change between qubits. Magnetic fields are used to bring specific qubits into 
resonance with the classical EM fields or the cavity mode. The two state qubits 
are coupled to a bath of EM field spontaneous emission (SE) modes, and the 
cavity mode is coupled to a bath of cavity decay modes. For large N the nu- 
merous vibrational modes of the qubits also act as a reservoir, coupled to the 
qubits, the cavity mode and the SE modes. 
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